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We calculate the optical conductivity of one-dimensional 
Mott insulators at low energies using a field theory descrip- 
tion. The square root singularity at the optical gap, char- 
acteristic of band insulators, is generally absent and appears 
only at the Luther-Emery point. We also show that only few 
particle processes contribute significantly to the optical con- 
ductivity over a wide range of frequencies and that the bare 
perturbative regime is recovered only at very large energies. 
We discuss possible applications of our results to quasi one- 
dimensional organic conductors. 

PACS numbers: 71.10.Pm, 72.80.Sk 



Measurements of dynamical properties and in particu- 
lar the optical conductivity cr(u>) are supposed to pro- 
vide a stringent test of the existing theories of quasi 
one-dimensional (ID) systems. The behaviour of cr(u) 
in the metallic regime is easily understood in terms of 
the Tomonaga-Luttinger theory 0] . The situation in the 
Mott insulating phase Q is much more complicated as 
a spectral gap is dynamically generated by interactions. 
Here u(u>) has until now only been studied by perturba- 
tive methods ffl5j, which are expected to work well at 
high and intermediate frequencies but are not applicable 
to the most interesting regime of frequencies close to the 
optical gap. The purpose of the present work is to deter- 
mine cr(w) in ID Mott insulators for all frequencies much 
smaller than the bandwidth, which is the large scale in 
the field theory approach to the problem. In particular 
we obtain for the first time the true behaviour of a(co) 
just above the optical gap. 

An important property of one-dimensional systems 
that significantly simplifies our analysis is spin-charge 
separation, which occurs at energies much smaller than 
the bandwidth. In this regime <t(uj) is determined solely 
by the charge degrees of freedom. The standard descrip- 
tion of the charge sector of ID Mott insulator is given by 
the sine-Gordon model (SGM) §,§ 



H sG = I dx 



47r(n) 2 



1 

167T 



(9^) 2 +2 Ai cos(/30) 



(1) 



Here the momentum and coordinate densities obey the 
standard commutation relation [n(a;), 4>{_y)\ — —i5{x—y). 
Throughout this letter we set the charge velocity and h 
equal to one. 

The cosine term in the Hamiltonian is related to Umk- 
lapp processes and the value of the sine-Gordon coupling 
constant (3 is determined by the interactions. The Umk- 
lapp processes are relevant for 1 < 1 and dynamically 



generate a spectral gap M, which is related to /j, by jlrj). 
For 1/2 < [3 2 < 1 the spectral gap is related to optical 
gap A (i.e. the gap seen in the optical absorption) by 
A = 2M whereas for /3 2 < 1/2 solitonic bound states are 
formed below 2M. 

Our calculations of cr(o->) are based on the exact solu- 
tion of the SGM and in particular on the work of Smirnov 
§. We confine our analysis to the repulsive regime 
1/2 < [3 2 < 1, where the excitation spectrum consists 
of charged particles and holes (solitons and anti-solitons), 
which do not form bound states. At the "Luther-Emery" 
point (3 2 — 1/2 the SGM is equivalent to the theory of 
free spinless massive Dirac fermions. In this limit the 
solitons become non-interacting particles and the Mott 
insulator turns into a conventional band insulator. In 
the limit 1 — > 1 the SGM acquires an SU(2) symmetry 
and describes the Hubbard model at half-filling in the 
regime of weak interactions |7|^] and ct(lj) was recently 
determined in 

The optical conductivity is related to the imaginary 
part of the current-current correlation function, x( w i <?) — 
(j-q 3q) b y 



<j(u > 0) = lm{x(uJ,q = 0)} /lo. 



(2) 



The current density operator is proportional to the mo- 
mentum density 



3q 



A^ 2 U n 



n„ 



dxll(t, x)e 



iqx 



(3) 



The non-universal coefficient A 1 / 2 depends on the de- 
tailed structure of the underlying microscopic lattice 
model. 

Using the spectral representation one can express the 
optical conductivity at T — as a sum over matrix ele- 
ments of the zero wave vector Fourier component of the 
momentum operator: 

*(w > 0) = - V \(0\IL \n)\ 2 S[uj - (E n - E )}. (4) 
iii — * 



Here |0) and |n) represent the ground state and excited 
states with energies Eq and E n respectively. The difficul- 
ties in computation of the optical response are related to 
the fact that one requires not only the knowledge of the 
spectrum E n , but also of the matrix elements of the mo- 
mentum operator. The exact expressions for the matrix 
elements (n|Ho|0) are extracted from the exact solution 
by means of the so-called form factor bootstrap procedure 
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Pj . This approach is particularly efficient for strongly in- 
teracting integrable models with spectral gaps, because 
for a given energy u> the spectral representation for the 
imaginary part contains only a finite number of terms 
(in the absence of bound states at most [o>/A] terms). In 
practice the spectral sum is found to converge extremely 
rapidly, so that a very good approximate description can 
be obtained by taking into account intermediate states 
with at most four particles The multi-particle ma- 
trix elements become essential only at very high energies 
where the field theory can no longer be used to describe 
the underlying lattice model anyway. 

In order to compute (^) we need to introduce a suit- 
able spectral representation. In the parameter regime 
we study, the spectrum contains only solitons and anti- 
solitons with relativistic dispersion e(j>) — \/p 2 + M 2 . 
It is useful to parametrize the spectrum in terms of a 
rapidity variable 9 



P 



M sinh / 



M cosh( 



(5) 



Solitons and anti-solitons are distinguished by the inter- 
nal index e = ±. A state of n solitons/anti-solitons with 
rapidities {9k} and internal indices {£k} is denoted by: 
\9 n . . . 0i) en ... si . Its total energy E, momentum P and 
electric charge Q are: 



P = M sinhflfe, E = M ^ cosh6> fc , Q cx ^e fe 



fc=i 



fc=i, 



fe=i 



In terms of this basis cr(u>) is expressed as 

*M = — L L J > m ■ ■ ■ 

n=0 Si J y ' 

x 6(M ^2 sinn 9k)$(u - M Y cosh ^fc) 

k k 
— <T 2 (uj) + (T4(cj) + ... 



Here 



Pi 



9n)ei...e 



(O|j(O,O)|0 n 



n/s n ... Sl 



(7) 



(8) 



are the form factors of the current operator, o<z (u>) and 
<ja{ui) represent the contributions from 2 and 4-particle 
processes and the dots indicate processes involving higher 
number of (anti)solitons. We note that as a consequence 
of symmetry properties only intermediate states with an 
even number of particles contribute to this correlation 
function. From (|7j) it is easy to see that only 2-particle 
processes contribute up to energies u> = AM, only 2 and 
4-particle processes up to w — 6M and so on. 

The form factors (^) have been determined in Q and 
can be used to calculate the first few terms in the expan- 
sion (0). We find 



„ 2M . ^ 9 'r?>wi 2 . 

lu 2 Vuj 2 - AM 2 

where Q{x) is the Heaviside function, 

2ttM sinh6>/2 
f(0) = P(0)+- = f(9)-+ = — 



(9) 



x exp ^ I dt 



ip C o S h(^) 

sinh 2 t(l - iO/n) sinh<(£ - 1) 
t sinh 2t cosh t sinh f £ 



(10) 



and 



9 = 2arccosh(cD) , £ = (3 2 /(l - 1 ) , u = w/2M. (11) 
The four particle contribution is of the form 



e(^4M) 



Ei <7=± 



f 6(9) 
d0 / d 7 

-6(9) 



x|/ J (5~ ^,S+^ff + + 7>S-» + 7) E i... E J S 




cosh 2 9 sinh 2 7 + a) 2 — cosh 6* cosh 7 I — 1 



c 2 r~ 2 



x [cosh 2 9 sinh 2 7 
where 

a = arccosh(fD — 1), b(9) = arccosh 
cosh(a/2) + exp (—7) cosh 9 



(12) 



2cosh6» 



(6) g = In 



a = 2arccosh 



cosh 9 sinh 7 + ui 2 — cosh 9 cosh 7 



The four particle form factor is given by 
4tt 3 

k<l 



/'(ft, ...,&)__++ 



3 £Md 2 nc(/3 fc -A) 

n (sinhp 2+m - /?„ - wr)^])" 1 



m,n=l ,2 

x2 sinh[(/? 4 +03 -Pi -fa- 2m)/2£] 



<exp(-i^,9 fe ) / ^J^v3(a-/? fc )cosh(a- i^/? fc ) 



x A (e 2Q /« I e 2ft /« , e 2& /« I e 2ft , e 2 ^/« ) . 



(13) 



The different orderings for (|13j) that appear in ([12]) can be 
obtained using the following property of the form factors 



f(0i 



■ ■ 9i+l 



+1 • ■ • 9 n ) El ... EiEi+1 ... en SJ J +1 {9i — 9 i+ i) 



(14) 



where S is the two body scattering matrix. The various 
functions appearing in Eqs ( |l3| ) and ( [l4| ) can be found in 
B (note that our definition of £ differs by a factor of tt) . 
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FIG. 1. Two particle (solind line) and one hundred times 
the four-particle (dashed line) contributions to the optical 
conductivity as a function of (u/M) for /3 2 = 0.9. Inset: 
threshold behaviour of the optical conductivity close to the 
Luther-Emergy point for four different values of /3; (3 = 0.75 
(solid), (3 = 0.74 (dotted), f3 = 0.73 (dashed) and fi = 0.72 
(long dashed) . 

The two and (one hundred times the) four-particle con- 
tributions to a(u>) for ft 2 = 0.9 are presented in Fig.l. 
Most importantly, the square root singularity, being a 
characteristic feature of band insulators, is suppressed 
by the momentum dependence of the soliton-antisoliton 
form factor and reappears only for the Luther-Emery 
point P 2 — 1/2. This effect was noted previously for 
the Hubbard model at half-filling |J which corresponds 
to the special SU(2)-symmetric point 1 — 1. We find 
that for any (3 2 ^ 1/2 there is a square root "shoulder" 
<t(uj) oc \Juj — A for u/A — 1 <C 1 as is shown in the in- 
set of Fig.l. In the vicinity of the Luther-Emery point 
/3 2 = 1/2 we obtain the following analytical expression 
valid for ui — 1 <C 1: 



<t{uj) oc — 



V^ 2 



1 



1] + £ 2 sin 2 7 



1 

2/32 



1 



(15) 



The square root singularity above u> = A for ft 2 — 1/2 is 
replaced by a maximum occurring at u>/ A — 1 oc j 2 . 

The four particle contribution to a is seen to be in- 
significant at low energies and becomes larger than the 
two particle contribution only at u> « 180M for (3 2 = 0.9. 
This suggests that the optical conductivity is well de- 
scribed by the combination of 2 and 4-particle contribu- 
tions up to several hundred times the mass gap. Com- 
putation of higher order terms in Eq.(j^) becomes cum- 
bersome and probably of no physical interest, since the 
previous analysis suggests that they become important 
outside the region of applicability of the field theory ap- 
proach to physical systems. 



At frequencies much larger than the gap it is possible 
to determine ct(uj) by perturbative methods. The leading 
asymptotics can be calculated by "conformal perturba- 
tion theory" jy]. Here the cosine interaction in (|l|) is 
considered as a (relevant) perturbation of the Gaussian 
model and correlation functions are calculated in a per- 
turbative expansion in powers of the scale /i, which then 
can be expressed in terms of the physical gap M as [n~2[ 



/' : 



r(/? 2 



7rr(l ~[3 2 ) 



M 



0Fr(l/2 + £/2) 
2r(£/2) 



(16) 



We find to leading order 



'<")= 28 - 4Sl (fwj) 2 "- ,lsI - 5> 



8tt 3 (3 2 



wr 2 (i-/? 2 )r 2 (i + /3 2 ) 



r(|) w 

2\/rf\±±«) M 



4/3-4 



(17) 



We emphasize that the ratio of the coefficients of the 
high- and low-energy asymptotics (0), (^|) is fixed ||, 
]l3[ . In other words, the amplitude of the power law 
in (|l7]) is tied to the overall factor in (^) and the form 
factor expansion must approach the perturbative result 
in the large-w limit. A comparison between the form 
factor results and (17) is shown in Fig. 2. We see that 



the asymptotic regime is not yet reached at energies as 
high as uj ~ 1000 M (in practical terms this implies that 
perturbation theory cannot be used to make contact with 
experiment). We note that the contributions due to inter- 
mediate states with 6,8,10 ... particles are all positive and 
will make the agreement of the form factor sum with per- 
turbation theory in the region to « 1000M only worse. A 
good way to overcome these deficiencies of bare perturba- 
tion theory is to carry out a renormalisation-group (RG) 
improvement as performed in Q. In Zamolodchikov's 
scheme |l^] the RG equations for the Sine-Gordon model 
are given by 



dg± 
dt 



9\\9x 



1 + 



The solution of (|T|) is 



<f<7|| 



.911 = 2 



9l 



l-(3 2 l + q 
P 2 1 - q ' 



where 



4(1 -P 2 ) 



(18) 



(19) 



(20) 



Using t- t = In (^57 
higher order terms as 



we can reexpress (|17|) up to 
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cr(w) = 



2wr 2 (2-/3 2 )r 2 (i + /? 2 ) 



r(|)e 3 / 4 ^ ' 

2 7/2 r (l+|) 



4/.r 



(21) 



The RG improved result (|T]) for cr(aj) is compared to 
the form factor result (sum of the two and four-particle 
contributions) in the inset of Fig. 2. The agreement is 
rather good down to energies of the order of 5M. 
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FIG. 2. Comparison between the 2 and 2+4-particle con- 
tribution to the optical conductivity and the perturbative re- 
sult, for f3 2 = 0.9. Inset: comparison between the form factor 
result and RG improved perturbation theory. 

One possible realisation of a ID Mott insulator are the 



(TMTSF) 2 X Bechgaard salts |14|. These materials are 
highly anisotropic and can be modelled as weakly cou- 
pled, quarter-filled chains. At energies or temperatures 
above the 1D-3D crossover scale E CI the interchain cou- 
pling becomes ineffective and a description in terms of a 
purely ID model with charge sector (Q) should be pos- 
sible ||. At present there is some uncertainty regarding 
the value of E cr because interactions can renormalize its 
bare value, set by the interchain coupling, downwards 
|jT5f . There is a lot of ambiguity in fitting our results 
to the data. The value of the optical gap 2M is not 
known and, as discussed above, we cannot calculate the 
overall normalisation of <j(uj). We therefore use these 
as parameters in order to obtain a good fit at large w 
(where the theory is expected to work best as 3D effects 
are unimportant) to the data fhij for any given value 
of 0. We obtain reasonable agreement with the data for 
(3 2 w 0.9, which corresponds to a Luttinger liquid param- 
eter of K p s» 0.23. This value is consistent with previous 
estimates (see the discussion in Q). 
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FIG. 3. Comparison between the optical conductivity 
calculated in the SGM for (3 2 — 0.9 (solid lines) and mea- 
sured optical conductivity for (TMTSF) 2 PF 6 from Ref.[12] 
(diamonds). The inset shows the same comparison on a log- 
arithmic scale. 

As is clear from Fig. 3, the model ([!]) seems to ap- 
ply well at high energies, but becomes inadequate at 
energies of the order of about 10 times the Mott gap 
(« 1600/cm in (TMTSF) 2 PF 6 ). Spectral weight is trans- 
ferred to lower energies and physics beyond that of a pure 
ID Mott insulator emerges. There are at least two mech- 
anisms that should be taken into account in this range 
of energies. Firstly, a small dimerization occurs in the 
ID chains and will almost certainly affect the structure 
of u{ijj) around its maximum. Secondly, the interchain 
hopping is no longer negligible [jl6| and ought to be taken 
into account. 

In summary, we have exactly calculated <j(ui) for a pure 
ID Mott insulator in a low-energy effective field theory 
approach. We have determined the threshold behaviour 
for the first time and found it to exhibit a universal square 
root increase for any (3 2 > 1/2. This is in contrast to the 
well-known suqare-root singularity that appears at the 
Luther-Emery point j3 2 = 1/2. In the "low" energy re- 
gion (w/A < 50) the optical conductivity is dominated 
by the two-particle form factor contribution with a small 
correction from four-particle processes. This means that 
the entire optical transport is dominated by two-particle 
processes! We furthermore have shown that the leading 
asymptotic behaviour obtained in perturbation theory is 
a good approximation only at extremely large frequen- 
cies, whereas RG-improved perturbation theory works 
well over a large region of energies. 
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